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1 Introduction 

The present experimental world average of the anomalous magnetic moment of the muon a^, 
assuming CPT-invariance, viz. = a^- , is 

^jexp) ^ 080 (63) X 10~" (0.54 ppm) , (1.1) 

where the total uncertainty includes a 0.46 ppm statistical uncertainty and a 0.28 ppm sys- 
tematic uncertainty, combined in quadrature. This result is largely dominated by the latest 
series of precise measurements carried out at the Brookhaven National Laboratory (BNL) by 
the E821 collaboration, with results reported in ref. [1] and references therein. The prediction 
of the Standard Model, as a result of contributions from many physicists is ^ 

^(SM) ^ ^4g^ ^ -^Q-ii ^ ^^ 2) 

where the error here is dominated at present by the lowest order hadronic vacuum polarization 
contribution uncertainty [4] (±42.0 x 10~^^), as well as by the contribution from the hadronic 

^For recent reviews see e.g. refs. [2, 3] and references therein. 
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light-by-light scattering, which is theoretically estimated to be (105 it 26) x 10^^^ [5]. The 
results quoted in (1.1) and (1.2) imply a significant 3.6 standard deviation between theory 
and experiment which deserves attention. In order to firmly attribute this discrepancy to new 
Physics, one would like to reduce the theoretical uncertainties as much as possible, parallel 
to the new experimental efforts towards an even more precise measurement of in the 
near future [6, 7]. It is therefore important to reexamine critically the various theoretical 
contributions to Eq. (1.2); in particular the hadronic contributions. Ideally, one would like 
to do that within the framework of Quantum Chromodynamics(QCD). Unfortunatley, this 
demands mastering QCD at all scales from short to long distances, something which is not 
under full analytic control at present. Therefore, one has to resort to experimental information 
whenever possible, to QCD inspired hadronic models, and to lattice QCD simulations which 
are as yet at an early stage. As a result, all the theoretical evaluations of the hadronic 
contributions to have systematic errors which are not easy to pin down rigorously. 

Our purpose here is to establish a simple reference model to evaluate the various hadronic 
contributions to within the same framework, and use it as a yardstick to compare with the 
more detailed evaluations in the literature. The reference model which we propose is based 
on the Constituent Chiral Quark Model (CxQM) [8]. This model emerged as an attempt to 
reconcile the successes of phenomenological quark models, like the De Riijula-Georgi-Glashow 
model [9], with QCD. The corresponding Lagrangian proposed by Manohar and Georgi (MG) 
is an effective field theory which incorporates the interactions of the low-lying pseudoscalar 
particles of the hadronic spectrum, the Nambu-Goldstone modes of the spontaneously broken 
chiral symmetry (S^SB), to lowest order in the chiral expansion [10] and in the presence of 
chirally rotated quark fields. Because of the SxSB, these quark fields appear to be massive. 
This model, in the presence of SU{3)l x SU{3)ii external sources has been reconsidered 
recently by one of us [11]. As emphasized by Weinberg [12], the corresponding effective 
Lagrangian is renormalizable in the Large-Nc limit; however, the number of the required 
counterterms depends crucially on the value of the coupling constant in the model and, 
as shown in [11], it is minimized for = 1. With this choice, and a value for the constituent 
quark mass fixed phenomenologically, the model reproduces rather well the values of several 
well known low energy constants. 

As discussed in ref. [11] the CxQM model has, however, its own limitations. Applications 
to the evaluation of low-energy observables involving the integration of Green's functions over 
a full range of euclidean momenta fail, in general, because there is no matching of the model 
to the QCD short-distance behaviour. There is, however, an exceptional class of low-energy 
observables for which the MG-Lagrangian predictions can be rather reliable. This is the 
case when the leading short-distance behaviour of the underlying Green's function of a given 
observable is governed by perturbative QCD. The decay — )• e"^e~, which was discussed in 
ref. [11], is one such example. Other interesting examples of this class of observables are the 
contributions to from Hadronic Vacuum Polarization, from the Hadronic Light-by-Light 
Scattering and from the Hadronic Zjj vertex ( provided, as we shall see, that the coupling g^ 
is fixed to gA = 1). The evaluation of these contributions with the CxQM Lagrangian is the 
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main purpose of this paper and they are discussed below in detail. They have the advantage 
of simplicity and can provide a consistency check with the more elaborated phenomenological 
approaches. 



2 Hadronic Vacuum Polarization. 

There is a well known representation [13] of the dominant contribution to the muon anomaly 
from the hadronic vacuum polarization shown in Fig. (1) 




Hadrons 

Figure 1. 

Hadronic Vacuum Polarization contribution to the Muon Anomaly. 



kg, - 2)hvp = a^") = r "^K (^) llmn(HVP)(i) , (2.1 



where ^ 



K^^\ = (A I i- ^,^\^,l^y (2.2) 



rn}, I vvr 



and ;^Imn(t) denotes the electromagnetic hadronic spectral function. It is a useful represen- 
tation because of the direct relation to the one-photon e+e~ annihilation cross-section into 
hadrons (rfig — )• 0): 

CT(t)|e+e-^(7)^hadrons} = ^ llmH^HVP) (t) , (2.3) 

and hence to experimental data, provided the necessary radiative corrections have been made 
to insure that one is using the one-photon cross-section. 



^The analytic expression of K (t/m^'j was first given in ref. [f4]; see also ref. [f9]. 
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In the CxQM with active u, d, s quarks and, to a first approximation, with neglect of 
gluonic corrections 



vr 



+ 



X 6 



6 



r 0{t - 4M^) , with 5 



AM? 



(2.4) 



The integral in Eq. (2.1) can then be easily done with the result shown in Fig. (3), the curve 
labeled (a), where the value for a|j^^^^ is plotted as a function of the only free parameter in 
the model, the constituent quark mass Mq ^ 

The constituent quark fields in the CxQM are assumed to have gluonic interactions as 
well but, since the Goldstone modes are already in the Lagrangian, the color-SU(3) coupling 
constant is supposed to be no longer running below a scale /Uq — 2 GeV where OsifJ-o) — 
0.33 and non-perturbative effects become significant. With inclusion of the leading gluonic 
corrections in perturbation theory, and to leading order in Large-Nc, the spectral function in 
Eq. (2.4) then becomes 



llmn(HVP)(t) 
vr 



vr/ 3 1 V2 6 
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N.asiVi) 3 



vr 
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PKsit) 



(2.5) 



where PKs{t) can be extracted from the early QED calculation of Kallen and Sabry [18] (see 
also ref. [19]): 
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(2.6) 



Also, at the level of the accuracy expected from the CxQM, it is sufficient to use the one loop 
expression 

asiVi) 1 



vr 



;^iv.-f)iog-^ 



with A ~ 250 MeV and nj = 3. 



(2.7) 



A 



There are many estimates of this contribution, as well as some of the higher order ones, with quark models 
which can be found in the literature. An earlier reference is [15] and two more recent ones [16] and [17]. 
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The resulting value for a/t in Eq. (2.1) in 10"^*^ units versus Mq in MeV is shown in 
Fig. (3), the curved labeled (b). 

In order to compare the CxQM results for a^^^^^ with the phenomenological determi- 
nations which incorporate experimental data, we still have to correct for the fact that the 
curve (b) in Fig. (3) only reflects the Large-Nc estimate of the model. As an estimate of 
the l/Nc^suppressed effects, we then consider the contributions from the 7r"'"7r~ and K~^K~ 
intermediate states to the spectral function in Eq. (2.1), as predicted by the CxQM. Notice 
that in this evaluation, the point like coupling {—ie){p^ — p'^) of scalar QED is replaced by 
the dressed coupling: 




Figure 2. 

Feynman diagrams contributing to the electromagnetic form factor T{Q'^) in Eq. (2.9). 



which, for = 1, is given by the expression: 



167r2 /2 I Q2 



dx X dy [1 — x(l — y)]- 



Ml 



Ml + xy{l - x)Q2 _ 



1 4Af2 



1 + ^ 



4M| 



1 + ^ + 1 



(2. 



In fact, this form factor, for gA 7^ 1, has UV-contributions which diverge and would require 
an explicit counterterm in the Lagrangian. The form factor in Eq. (2.9) has the asymptotic 
behaviour: 



lim F{Q' 



levrV- 



3 + 30 M2+^1^M^ 



(2.10) 
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and, in particular, fixes the value of the coupling constant Lg in the xPT effective Lagrangian 
to [20]: 

Lg = ^-. (2.11) 

Also, for = —t and t > 4Mq the form factor develops an imaginary part: 




-4M2), (2.12) 
and the form factor F{Q'^) obeys a once subtracted dispersion relation 

HQ^) = r , . - li^m , (2.13) 

J AMI 



t + Q^ -ie t J TT 



the subtraction ensuring that J-{Q'^ = 0) = , as fixed by lowest order xPT. 

The form factor J^{Q'^), however, does not match the QCD behaviour at large-Q^ values 
and, therefore, the estimate we propose for the 1 /Nc-suppresed contributions to the the muon 
anomaly can only be considered as reasonable up to values of t in Eq. (2.1) below t ~ /iq 
where the asymptotic pQCD regime sets in. Contributions beyond t ^ fiQ have already been 
taken into account by the second term of the spectral function in Eq. (2.5). 

The total contribution to a|i^^^^ in the CxQM, which incorporates gluonic contributions 
in the spectral function in Eq. (2.5) as well as the subleading Tr~^TT~ and K~^K~ contributions 
in the way described above is shown in Fig. (3) as a function of Mq, the curve labeled (c). 

These considerations provide us with a framework to fix the constituent quark mass Mq. 
The prediction of the CxQM, as described above, should be compared to the phenomenolog- 
ical contribution from hadrons formed of u, d and s quarks only, at the level of one-photon 
exchange. Contributions like for example the one from an intermediate state should there- 
fore be excluded so far (more on that later on), as well as those involving c, b and t quarks. 
From the numbers quoted in TABLE II of ref. [4], we then find that this restriction reduces 
the phenomenological determination of the anomaly from hadronic vacuum polarization to a 
central value 

4"^P)|phcn. ^ 653 x 10-10 (2.14) 
which, when compared with the results plotted in Fig. (3), shows that fixing Mq in the range 

Mq = (240 ± 10) MeV , (2.15) 

reproduces the phenomenological determination within an error of less than 10%. This deter- 
mination of the constituent quark mass is the value which we shall systematically use for Mq 
when evaluating the predictions for the other hadronic contributions to the muon anomaly. 
We shall then compare them to the various phenomenological determinations in the literature. 

We wish to emphasize, however, that the error of 10 MeV in Eq. (2.15) only reflects 
the phenomenological choice that we have made in order to fix Mq. As discussed in the 
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Mq (MeV) 
Figure 3. 

^(HVP) CxQM. Curve (a) is the contribution using the spectral function in Eq. (2.4); curve 

(h) the contribution using the corrected spectral function in Eq. (2.5) and curve (c) the contribution 
using the corrected spectral function in Eq. (2.5) with subleading tt+tt^ and K^K^ contributions 
incorporated as discussed in the text. 



Introduction the CxQM is only a model of low energy QCD and, as such, there is no a priori 
way to fix Mq from first principles. The error in Eq. (2.15) does not reflect the systematic 
error due to other plausible ways of fixing Mq. 

At this stage we wish to point out that the recent lattice QCD determination of a^f^^^^ 
with two flavours reported in ref. [21] can also be very well digested with a value of Mq within 
the range given in Eq. (2.15). 



3 Hadronic Vacuum Polarization Contributions at Next to Leading Order. 

The Hadronic Vacuum Polarization contributions at O (^)^ were classified long time ago in 
ref. [15]. Let us discuss their evaluation in the CxQM. 

3.1 Class A: HVP INSERTIONS IN THE FOURTH ORDER QED VERTEX DIAGRAMS. 

They correspond to the Feynman diagrams shown in Fig. (4), where the diagrams in each line 
in this figure are well-defined gauge invariant subsets. Here, the equivalent of the function 
K {t/m?^ in Eq. (2.2) at the two loop level, which we call iiT^'^^ (t/m^), was calculated ana- 
lytically by Barbieri and Remiddi [22]. The exact expression is, however, rather cumbersome 
and for our purposes it is more convenient to use an expansion of this function in powers of 
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-p, which is justified by the fact that the hadronic threshold in the integral that gives the 
contribution from the diagrams of Class A: 

a(HVP-A) = r ^K(^)(t/ml) llmn(HVP)(i) , (3.1) 

starts at 4:Mq ^ m^. The terms in the expansion in question for the kernel K^^^ {t/nifj) 
which we have retained are: 




Using the CxQM spectral function in Eq. (2.5), we find for this contribution the following 
result: 

^(HVP-A)^^^ = 240 MeV) = -171 x 10"" , (3.3) 

with a range 

- 181 X 10"^^ < aj,"^^"^) < -161 X 10"" , for 230 MeV < Mq < 250 MeV . (3.4) 

This result is to be compared with the phenomenological determination [23]: 

^(HVP-A) ^ _(207.3 ± 1.9) X 10-" . (3.5) 

We conclude that the CxQM reproduces, within the expected accuracy of the model, this 
phenomenological value, specially if we take into account that the phenomenological determi- 
nation includes contributions subleading in 1/Nc and from higher flavours, which are beyond 
the duality domain of the model. 

3.2 Class B: HVP insertions in the QED vertex with an electron loop. 

This is the contribution from the two Feynman diagrams in Fig. (5). A convenient represen- 
tation [15] for this contribution, is the one given by the integral 
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Figure 5. 

Feynman diagrams corresponding to the Class B contribution to the muon anomaly. 

denotes the real part of the electron self-energy. Using the CxQM spectral function in 
Eq. (2.5), we find for this contribution the following results: 

^(HVP-B)^^^ = 240 MeV) = 88.9 x 10"^^ , (3.9) 

with a range 

82.6 X 10"^^ < aj,^^'""^) < 95.9 X 10"" , for 250 MeV > Mq > 230 MeV . (3.10) 

This result is to be compared with the phenomenological determination of this contribu- 
tion which gives [23]: 

^(HVP-B) ^ (^Qg Q ^ Q ^ _ ^3 -^^^ 

We find again that, within the expected accuracy of the model, the CxQM reproduces the 
phenomenological determination. 

3.3 Class C: Iterated HVP Contributions. 



X 




Figure 6. 

HVP Contributions at 0{ol). 
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This is the contribution in Fig. (6) induced by the quadratic term in the expansion of the 
photon propagator in the lowest order vertex, fully dressed by hadronic vacuum polarization 
corrections, i.e. 

Do.p = [9., -^)-2 ^^ n(HVP)(g2) - ' (3-12) 

where Il^^^^\q^) denotes the proper vacuum polarization self-energy contribution induced 
by hadrons and a is a parameter reflecting the gauge freedom in the free-field propagator 
(a = 1 in the Feynman gauge). In fact, since the diagrams we are considering are gauge 
independent, terms proportional to qaqp do not contribute to their evaluation. The lowest 
order muon anomaly is then modified as follows: 

-\-^(-\[\x{l-x) ^, (3.13) 

n)2 KnlJ, \ + n(HVP) (^^2) ' 

and the perturbation theory expansion generates a series in powers of the self-energy function 



i + n{HVP) (z 



1 _ n(HVP) .z^ml ] + 

x± 2] \1 X 

-X M 



2 



n(HVP) -^ml 
l-x ^ 



2 



+ •••. (3.14) 



Writing a dispersion relation for each power of n^^^^-* to lowest order in the elec- 

tromagnetic hadronic interaction, i.e., 



n(HVP) {^mA = r ^ ^"f . ilmn(HVP) (t) , (3.15) 

results then, from the quadratic term of the expansion in Eq. (3.14), in the following repre- 
sentation [15] for the contribution from the Feynman diagrams in Fig. 6 



^(HVP-C) 



r ^ ilmnCVP) (t) r ^:i„„n<"VP) (i) k[±,.JL\ (3.16) 



with 

K{^,^] = i-) dx-. ^' 3.17 

\mf^ mf^J KttJ Jo U2 + A^^i - x)] (x^ + ^{1 - x) 



a composite kernel which correlates the two spectral functions. Using the CxQM spectral 
function in Eq. (2.5) for both ilmn(HVP) (t) and ilmn(HVP) (^f), we find a small contribution 
from this C-class: 

aJf^P-^)(MQ = 240 MeV) = 2.2 x 10"^^ , (3.18) 

with a range 

1.9 X 10"^^ < aj^^^P-^) < 2.5 X 10"" , for 250 MeV > Mq > 230 MeV . (3.19) 
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Two independent phenomenological determinations of this contribution (with errors which 
are likely to have been underestimated) give: 



^(HVP-C)(j.gf^ [23]) = (3.4 ±0.1) X 10' 



11 



and 



amP-C)(rei. [30]) = (3.0 ± 0.1) x lO'^^ 
Again, they compare reasonably well with the CxQM prediction. 

• Why is this contribution so small? 



(3.20) 
(3.21) 



This is an interesting question which, to our knowledge, has not been addressed in the 
literature. We wish to take the opportunity to answer it here. 

The main point is the following: instead of writing a dispersion relation for each power of 

n(HVP) (^2) 

, we could have chosen to write a dispersion relation for the squared photon 
self-energy [n(^^^)(g2)] , i.e. 



n(HVP)(^2 



dt 



t t 



2Ren(HVP)(t)iiinn(^iVP)(t). 



This leads to a representation for the muon anomaly {q^ = jzi^rn?^): 



rnvp-c) ^ ^ 



dx 



a;2 + ^(l-x) 



-2Ren(HVP)(t) -Imn(HVP)(t) 
vr 



(3.22) 



(3.23) 



similar to the one we have used in Eq. (3.6) for the evaluation of afi^^ Gauge 
invariance guarantees that the subtraction constant in the double dispersion relation 



r^-^- ilmn(HVP)(t) X ^\ . ilmn(HVP)(t') (3.24) 

Jo t t — q^ — le TT Jo r r — g — ze TT 

and the one in the single dispersion relation in Eq. (3.22) are the same, so that the 
physical electric charge corresponds to the one measured classically. In other words, 
gauge invariance guarantees that the physical content of the two equations (3.24) and 
(3.22) must be the same. Yet, algebraically, starting with the r.h.s. in Eq. (3.24) and 
using the partial fraction decomposition: 

1 1 1111 , , 



t — q^ — ie t' — q^ — ie t — q^ — ie t' — t t' — q^ — ie t — t' 
one gets the following relation: 

n(HVP)(<^2) ^ n(HVP)(g2)= fn(HVP)(g2)l'_g2 2Ren(HVP)(i)limn("vp)(i) . (3.26) 

L J Jo * TT 
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Obviously, the only way to preserve the identity n^^^^) {q'^) x n^^^^) (q^) = [n^^^^) (q^)] ^ 
is that 

r % 2 Ren(HVP)(t)ilmn(^VP)(^^ ^ q ^ (3 27) 

which is a highly non trivial constraint! ^. It is this constraint which answers the ques- 
tion of why Oy^^^ turns out to be so small. Indeed, it implies that the a priori 
leading term of 0{m?^) in an expansion in powers of in the r.h.s. of Eq. (3.23), con- 
trary to what happens with the lowest order hadronic vacuum polarization contribution 
in Eq. (2.1) where it provides the dominant contribution, is not there in the dou- 
ble hadronic vacuum polarization contribution. The leading term in a m^-expansion 

for a^^^^ must be 0{m't) at least. In fact, a detailed analysis shows that it is 



log 



with Mil a hadronic scale which in the CxQM is Mq of course. This 
is the reason why the double hadronic vacuum polarization contribution to the muon 
anomaly is so small and, as we have shown, this is a model independent statement. 

Comment on Radiative Corrections 

Hadronic vacuum polarization generates part of the radiative corrections to the total 
e~^e~ annihilation bare cross-section into hadrons. In fact this correction leads to the 
following modification of the bare cross-section 



CT(t)e+e-^hadrons = ^ -lmn(HVP) (t) 

t vr 

2 Ren(«VP)(t)l iimn(HVP) (t) , (3.28) 
t I J vr 



which corresponds to the modification 

ilmn(HVP) (t) ^ h _ 2 Ren(HVP)(t)l ilmn(HVP) (t) , (3.29) 
vr L J vr 

and leads, precisely, to the muon anomaly contribution given in Eq. (3.23). In other 
words, if in the lowest order expression for the muon anomaly one inserts the bare total 
e~^e~ annihilation cross-section into hadrons, we are indeed calculating the lowest order 
contribution o|j^^^^ in Eq. (2.1). This implies that the appropriate radiative corrections 
to the physical cross-section have been made including the correction due to hadronic 
vacuum polarization. The alternative is to leave the physical cross-section uncorrected 
for hadronic vacuum polarization, in which case, when inserted in the lowest order 
expression, one is then calculating: a^^^^^ -|- a|j^^^ . Then, obviously, one should not 
add an extra independent evaluation of a^^^^ . 



In QED, one can easily check this constraint in perturbation theory. 



- 13 - 



The warning here, specially for theorists, is that in using experimental hadronic cross- 
sections to compute hadronic vacuum polarization contributions to the muon anomaly, 
one should be very careful to know exactly what these cross-sections correspond to. 
Often the data which is used corresponds to different experiments which complicates 
even further the issue. 

Another warning, this one for experimental physicists, concerns the dynamical con- 
straint given in Eq. (3.27). In doing hadronic vacuum polarization corrections to the 
total cross-section numerically (e.g. involving an iterative procedure, as mentioned in 
some of the experimental papers) one should be careful to check that this constraint, 
which involves rather subtle cancellations, is indeed satisfied. 

3.4 Class D: Contributions with HVP corrections at 0{a) 



X 




Figure 7. 

Feynman diagrams corresponding to the Class D contribution to the muon anomaly. 

In the CxQM there are two types of contributions to this class: the n^'y exchange and the 
constituent quark loop with a virtual photon insertion. They correspond to the photon 
propagator content illustrated in Fig. (8). Corresponding to these two subclasses we shall 
write 

and discuss separately the two contributions. They are both leading in the Large-Nc limit. 
3.4.1 Contribution from the tt'^j intermediate state 

Here it is convenient to use the representation (see page 231 in ref. [24] and ref. [25] ^): 




(3.31) 



where Il^'^°'^\k'^) denotes the renormalized photon self-energy from the contribution and 
the integration is over the value of the self-energy in the euclidean. In fact, we find that a 

^This is a representation which is now often used by our lattice QCD coUeagues to compute the hadronic 
vacuum polarization contribution. 
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vwM^prwvA/ 




+ 



(b) 



Figure 8. 

Feynman Diagrams which contribute to the Photon Self-Energy to 0(0^) in the CxQM. 



better representation, which avoids renormahzation issues, is the one in terms of the Adler 
function [26] ^. Using integration by parts in Eq. (3.31) with (1 — z) = —^^(1 — z)"^ and the 
fact that n(0) = 0, one finds 



with the Adler function (Q^ = -k'^) 



-m, 



(3.32) 



^(-°7)(Q2)^_g2 



dn(^°^)(A:^) 



12 _ ^ 2 

1-z ^ 



(3.33) 



In the CxQM, the ■K^'yy three-point function at each vertex in the first diagram of 
Fig. 10(a) can be expressed in terms of the following parametric representation: 



-le 



/ dxx 
Jo Jo 



dy 



12^2/. 
2M^ 



Mq — x(l — x)(l — y)k'^ — x'^y{l — y){k + qY — xy{l — x)q'^ — ie 



(3.34) 



Here, the constituent quark mass Mq acts as an UV-regulator of the contribution to the 
muon anomaly. In the limit Mq — t- oo this form factor reduces to the 7r'^77 Adler, Bell-Jackiw 



Fred Jegerlehner has often advocated also the advantages of this representation. 
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point-like coupling (ABJ): 



A/o— S'oo 



2 



(3.35) 



and in this limit, the contribution to the muon anomaly becomes UV-divergent. Using 
dimensional regularization and the MS-renormalization scheme, the result in this limit, with 
<C 771^ for further simplification, and ji the renormalization scale, is 



log^ + § + 
mt 



-11 



mt 



In particular, for /i = Mp, one finds 

a(:"^)(A5J)p=M, = 2.5x 10 
a result which, within a 30% error, is consistent with the one in ref. [27]: 

a^^^^ ~ 3.7 X 10-" , 



m 



obtained with Vector Meson Dominance like form factors: 



(3.36) 

(3.37) 
(3.38) 

(3.39) 



127r2/^ M2 - A;2 M2 - g2 • 

The Feynman parameterization of the full tt^^ Adler function in the CxQM, using the 
form factor expression in Eq. (3.34), results in the following representation: 



^ N?: — ^r-nr - [ duv? [ dvv [ dw 

.TT/ ^167r2/2 9 io io io 

X / dxx j dy I dx'x' [ dy' — — 

/o Jo Jo Jo xy(l - xy)x'y'{l - x'y') 



y{l - xy) 



y'{l - x'y') 



+uv{l — w) — [uv{l — w) + u{l — v)6 + (1 — u)(5'] ^1 



7e2 



(3.40) 



where 



^2 ^2 f(i-y)[i-a;(i-y)] . n {'^-y')['^-x'(i-y')], 

n^ = Q^ A 1 -v) + ^ ; ''^ {i - u) 

[ y(l - xy) y'[l - x'y') 

+uv{l — w) — \uv{l — w) + u(l — v)5 + (1 — n)(5'] ^1 



+ 



l-u 



u{l — v) 
xy{l — xy) ' x'y'{l — x'y') 



+ m\uv{l — w) . 



(3.41) 
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and 

S= '^'-y^ and 6' = -f} . (3.42) 
1 — xy 1 — x'y' 

Performing the integration over the seven Feynman parameters numerically we find the fol- 
lowing results: 

aj7^) {Mq = 240 MeV) = 2.17 x 10~" , (3.43) 

with a range 

2.10 X 10~^^ < aj;'^) < 2.18 X IQ-^S for 230 MeV < Mq < 250 MeV , (3.44) 

consistent with the estimate in Eq. (3.37). 

We observe, however, that these results turn out to be an order of magnitude smaller 
than the phenomenological contribution quoted in the TABLE II of ref . [4] (see also ref . [23] ) 

a|7T) = (44.2 ± 1.9) X 10"" , (3.45) 

which uses as input the measured a{e'^e~ — ?• n^'^) cross section in the energy interval 0.60 < 
^ < 1.03 GeV [28]. 

• Why this discrepancy? 

In order to understand better the underlying physics let us use instead the representation 
for aj^'^'' in terms of the spectral function ^\mn'^'^^\t) i.e., 



r -K ( A,] ilmn(-T) (t) . (3.46) 



The phenomenological determinations of af^'~^^ in refs. [4, 23] implicitly assume that 
-lmn^'^'^^(t) is completely saturated by the intermediate state. Notice however that 
in the CxQM there are other intermediate states which also contribute to ^Imn('^''')(t); 
they correspond to the QQ-, QQ^ and QQtt discontinuities in the diagram (a) of Fig. (8). 
These discontinuities are automatically included in the calculation of a!j^ '^'^ which uses 
the euclidean representation in Eq. (3.32). In order to compare the CxQM determina- 
tion to the phenomenological ones, let us then restrict ;^Imn^'^''') {t) to the contribution 
from the on-shell intermediate state only. Then 

llmn(--)(t)| = t fl - -F^^^r^ (m^,t,0) " , (3.47) 

with 



J^^^^,^) {ml,t, 0) = / ^ dxx C dy^ '^^^ , , . (3.48) 
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This form factor has an imaginary part which, for m-p: < 2Mq, is: 

4M2 



and a real part: 



t-ml 2M, 



log 



i + \n 



t 



eit-iM^). (3.49) 



t — ml 



arctan 



Q 1 



+ < 



arctan 



In' 



2 

2A/o 



1 + VI 



ml<t< 4M^ 



(3.50) 



The shape of the spectral function in Eq. (3.47), in units of (^)^ and for the value 
Mq = 240 MeV, is shown in Fig. (9). 



1 



-Imn"'''(t) 0.10 



t (GeV') 



Figure 9. 

The Spectral Function in Eq. (3.47), in units of (f )^ for Mq = 240 MeV. 



The result of the integral in Eq. (3.47) with the spectral function plotted in Fig. 9 is 
then 

a{7^)(Eq. (3.46)) = 3.0 x 10"^^ for Mq = 240 MeV , (3.51) 
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a result which is shghtly higher than the full CxQM contribution in Eq. (3.43) but still 
well below the phenomenological determinations [4, 23]. 

Let us then try to simplify the phenomenological determination as much as possible to 
see where the big contribution comes from. For that, it will be sufficient to approximate 
Eq. (3.46) as follows: 



2 

(7r7) 



vr/ r o t vr 



and use a narrow width expression for the spectral function, which as we shall soon see, 
is dominated by the oj contribution. This results in the simple formula: 

3M2 7r M,, r,, ' ^ ' 



which reproduces, in order of magnitude, the phenomenological estimates. We can, 
therefore, see that the hig number comes from the large experimental value of the 
branching ratio 

^ ~ 8 X 10~2 . (3.54) 

t UJ 

Notice that in the case of the p contribution the corresponding branching ratio is much 
smaller: 

Tip 7r°7) 4 

^ ~ 6 X 10"^ . (3.55) 

p 

It is the large branching ratio in Eq. (3.54) which the C^QM fails to reproduce! 

Phenomenologically, the large branching ratio p^'^Z^^o^j is due to the a;-0 mixing and 
the fact that the (\) is an almost pure ss state By construction, the CxQM form 
factor is SU{3) invariant and, therefore, like any model which is SU{3) invariant, fails 
to reproduce this phenomenological fact. 

We are aware of the fact that in the CxQM there are also further contributions of the 
71*^7 subclass: those from the r/87 and 7707 intermediate states. We refrain from discussing 
them because their comparison with their corresponding phenomenological determinations 
requires issues like 773 ~ Vo mixing as well as the question of the rj' mass in Large-Nc which 
are beyond the scope of the model we are discussing. 

3.4.2 Contribution from the Quark Loop to 0{a) 

This contribution is given by the following integral representation: 

= / ^K{t/ml) -Imn(« it) , (3.56) 



^We thank Marc Knecht for reminding us of this. 
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where 



llmnW'")(t)=f- 



TT 



+ 



+ 



PKs{t) , 



(3.57) 



with pKsit) given in Eq. (2.6). As Mq — t- oo it decouples with a leading behaviour: 



9 3M2 



vr 



dtM^ 



41 



, , PKs{t) = — 
4M2 t t 162^ 



13.7 X 10~" for Mq = 240 MeV. 



(3.58) 



The full numerical evaluation gives 



a 



{Mq = 240 MeV) = 11.4 x 10"" , (3.59) 



with a range 

10.6 X 10"" < ajf?'") < 12.3 x 10"" , for 250 MeV > Mq > 230 MeV . (3.60) 

Altogether we find that, although the exchange contribution increases logarithmically 
as a function of Mq, while the quark loop decouples as an inverse power of Mq, their ratio 
for Mq large goes as 

M^ 27 
Mq~^oo Wtt^'j^ 82 



and, therefore, for values of the constituent quark mass in the range 250 MeV > Mq > 
230 MeV, it is the quark loop contribution which still dominates. 

The total sum of the two contributions of Class D in the CxQM is then: 

^(HVP-D)^^^^ = 240 MeV) = 13.6 x 10"" , (3.62) 

with a range for this total 

12.8 X 10"" < aj,^^^"^) < 14.4 X 10"" , for 250 MeV > Afg > 230 MeV . (3.63) 

Except for the contribution, it is difficult to compare the overall CxQM prediction for 
the Class D contributions with the phenomenological estimates. The reason is that, a 
priori, when inserting a physical observable to evaluate the diagram in Fig. (7) one needs 
two types of contributions: the one from the cross section a{e~^e~ — t- Hadrons + 7) and 
the one from the interference of the amplitude e~^e~ — )• Hadrons with the same amplitude 
where a virtual photon has been emitted and reabsorbed. In fact, individually, these two 
contributions are infrared divergent, which complicates things even more. This is a place 
where it would be interesting to see if lattice QCD can eventually make an estimate of these 
Class D contributions which, so far, remain poorly known phenomenologically. 



- 20 - 



Table 1 gives a summary of the results for the four classes of contributions discussed here 
and evaluated within the framework of the CxQM, with a total 

^[HVP-(A,B,c,D)] ^ (_g4 ^ X 10^1^ . (3.64) 

This CxQM result is to be compared with the number quoted in the latest evaluation 
in ref. [29] for this contribution which, however, does not include the important contribu- 
tions from the and ?/7 intermediate states already incorporated in the lowest order HVP 
contribution: 

^[HVP-ncxtordcr]^g+g-) = (-98.4 ± 0.6cxp ± OA,^) X 10"^^ . (3.65) 

Again, except for the n^'y issue already discussed, the agreement within the errors of the 
model is quite reasonable. 

Table 1. Results for the HVP contributions of O {^f in the CxQM. 



Class 


Result in 10 units 


A 


-171 ± 10 


B 


89 ± 7 


C 


2.2 ±0.3 


D(7r7) 


2.2 ±0.1 


D(Q-loop) 


13.5 ±0.5 


Total 


-64± 12 



4 Hadronic Light by Light Scattering Contributions. 

The standard representation of the contribution to the muon anomaly from the hadronic 
light-by-light scattering shown in Fig. (10) is given by the integral [31]: 

^ tt(H) / ^" 

Jg=0 

xtr j(^±m^)[7^7^](^+m^)7^-— i Y-T-T T^'') ' (^'l) 

where Uliupcriq, qi,q3, q2), with q = p2—pi = —qi—q2—q3, denotes the off-shell photon-photon 
scattering amplitude induced by hadrons, 



(HLbyL) ^ -^e / d-^qi f ^ g2 1 

48m. J (27r)V (27r)4 q^q^q, + q^V 
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Figure 10. 

Hadronic Light-by-Light Scattering contribution to the Muon Anomaly. 



X (0|r{ J^(0) , J,(xi) , Jp(x2) , J.(X3)}|0) , (4.2) 

and J^{x) = QqQ{x)^nq{x) is the Standard Model electromagnetic hadronic current where, 
for the light quarks, Qq = diag(2/3, —1/3, —1/3). 

In the CxQM there are two types of contributions: the Constituent Quark Loop (CQL) 
contribution shown in Fig. (11) and the Goldstone Exchange Contribution shown in Fig. (12) 

X 




Perm. 



Figure 11. 

Constituent Quark Loop Contribution to the Muon Anomaly in the CxQM. 



with constituent quark loops at each vertex. We shall consider these two types of contribu- 
tions, both leading in the 1/Nc-expansion, separately. . 

4.1 Class A: The Constituent Quark Loop Contribution. 

This contribution can be obtained from the QED analytic calculation of Laporta and Rem- 
midi [32] with the result 
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Figure 12. 

Goldstone Exchange Contribution to the Muon Anomaly in the CxQM. 



a 



(HLbyL) 



(CQL) 




2^^ ^ 16 



m. 



161 , 2 ^^O 16189 , Ml 13 161 vr^ 831931 

log -ZT + tt;;:^ log — r " — C(3) + ttttttt— + 



3240 



97200 



mf 



18 



1620 6 972000 



+ O 



m 



m. 



log 



(4.3) 



A plot of this contribution versus Mq is shown in Fig. (13). We find 

^(HLbyL) ^QQL) ^ 82.2 X 10~" at Mq = 240 MeV , (4.4) 

with a range 

76.3 X 10"^^ < aj.^^'^y^) (CQL) < 88.8 x 10"^^ for 250 MeV > Mq > 230 MeV . (4.5) 

The gluonic corrections of O (^) to the leading term in Eq. (4.3) have been recently 
calculated in ref. [33] and found to be rather small. 

4.2 Class B: The 7r° Exchange Contribution. 

The expression for this contribution in terms of the vertex form factors F^o*^*^* and F^o*^^* 
can be found in ref. [34]. When applied to the CxQM we have: 



„(HLbyL)/ On = p2 ^ f d'^Ql [ dS}2 

JxQM 3 ] (2^)4 ] (2^)4 g2g2(^^ + ^^)2[^2 + 2p . q,][ql - 2p • q^] 



^1(91,92;^)- 



Q2 - "^l 



+ T2{qi,q2;p) 



•^ff,^^, ((gi + g2)^ ql gj) 4ff;f; ((gi + g2)^ (gi + g2)^ 0) 



(gi + 92)^ - m2 



(4-6) 
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where 



Ti{qi,q2;p) = 'i{p ■ qi){p ■ q2){qi ■ 92) - 2(p • q2fql - (p • qi){qi ■ ^2)^2 

+ 3(p • q2)qlql - 2{p ■ q2){qi ■ q2? + 2ml [ll^l " (^l " ^2)'] (4.7) 

originates in the first and second diagrams of Fig. (12), which give identical contributions, 
while 



T2{qi,q2;p) = 2{p ■ qi){p ■ q2){qi ■ ^2) - 2{p ■ qifql + {p ■ qi){qi ■ 92)^2 

+ {p ■ qi)qlql + ^ [^ihl - (qi ■ Q2f] (4.8) 

originates in the third diagram of Fig. (12). 

It is well known [35] that, asymptotically for Mq ^ m^r, the vr'^-exchange contribution 
must behave as: 

»r-'(.\«M ^ OV^i£%[^"=^-o(>"^) . (4..) 

In the CxQM this contribution can be evaluated exactly. Notice that J'^?^}^^ {q^,q'^,0) has 
a simple analytic expression: 
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-^e 



127r2/, 



dx- 



Ml 



le 



1 + ^ 



:l0g 



1 + ^-1 



1 + ^ + 1 



-le 



.2 Ne 1 

127r2/7r27ri 7^ 



(is 



[-Ml] i r(2 - 2s) 



(4.10) 



where the expression in the third hne gives a very useful representation for analytic evalua- 
tions. The full expression of the other vertex function is: 



Tr(xQM) / 2 2 2\ _ -2 

1 /■! 2M2 

'^'^'^Jo '^^M^ - x{l -x){l- y)ql - x^y{l - y)ql - xy{l - x)ql ' ^^'^^^ 



for which the following representation can also be used: 



Tr(xQM) / 2 2 2\ _ -2 

ci— ioo C2—ioa C3— ioo 

r(l - Si - S2)r(l - Sl - S3)r(l - S2 - Ss) 



r(3-2si -2s2 -2s3) 



r(si)r(s2)r(s3)r(i - si - S2 - S3) • (4.12) 



The interest of this Mellin-Barnes representation is that it only modifies the powers of the 
propagators: q\, ™d {qi + (72)^ in the first line of Eq. (4.6), and it provides a systematic 

2 2 

way ® to compute the asymptotic expansion in ^Jr ^-nd —3- powers, and powers of logarithms. 
We postpone, however, this analytic calculation to a forthcoming publication and, instead, 
proceed here to a numerical evaluation. 

In order to evaluate a^^^^^ {t^^)xCim in Eq. (4.6) numerically, it is useful to apply to the 
integrand in that equation the technique of Gegenbauer polynomial expansion, as was done 
in ref. [34]. Then one can reduce the q\ and §2 integrations to two euclidean integrals over 
Q\ = —q\ and Q\ = — ( both from to 00), and an integral over cos 6 with 6 the angle 

*See e.g. ref. [36] for an example of this method. 
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between the two euclidean four-vectors Qi and Q2- The integrand in question, which is 
exphcitly given in ref. [3], is then very convenient for numerical integration. 
We find 

mLbyL)^^o^^^^ ^ gg Q ^ ^Q-ii fQj. Mq = 240 MeY, (4.13) 



with a range 



and 



The result 



5,^^^y^)(7r°)^QM = 64.6 X 10"^^ for MQ = 230 MeV, (4.14) 
aj,^^''y^)(^°);^QM = 71.3 X 10"" for AfQ = 250 MeV. (4.15) 
^(HLbyL)(^0)^^^ = (68 ± 3) X 10-11 , (4.16) 



which does not include the systematic error of the model, agrees well with the phenomenolog- 
ical determinations of this contribution which, according to the most recent update [38] and 
depending on the underlying phenomenological model for the form factors J-^?^]^^ [q2, Qi^Qs) 
vary between 

^(HLbyL)^^O)^^^^ = (57.4 ± 4.6) X 10-11 and a'-^^^^^\-K\^^^, = (80.1 ± 4.7) x lO-n . 

(4.17) 

Again, for the same reasons mentioned at the end of Section III.4.1, we do not discuss 
here the contributions from the rj and t]' exchanges. 

It is a fact that asymptotically, for Mq — )• oo, the vr'^ contribution largely dominates the 
Constituent Quark Loop contribution: 

,{HLbyL)^^0^ N,M^ 1 . 



'(vrU) 



Mg^oo 167r2/2 C(3) - i rn. 



log' ; (4.18) 



however, this asymptotic behaviour is far from being reached at values of Mq between 
230 MeV and 250 MeV, for which the Constituent Quark Loop contribution still dominates 
over the Goldstone contribution. 

For the total hadronic light-by-light contribution in the CxQM, which includes the quark 
loop contribution as well as the vr'^-exchange contributions, we then find 

148 X 10-11 < a|,iii^''yi^)(CxQM) < 153 x lO-n , for 250 MeV > Mq > 230 MeV . (4.19) 

This result, which does not include the systematic error of the model, has to be compared 
with the phenomenological estimate 

^(HLbyL) ^ ^^22 ± 18) X 10-11 , (4.20) 

for the total of the hadronic contributions not suppressed in the 1/Nc-expansion (see e.g. 
ref. [39] for details). Within the expected systematic uncertainties they compare rather well. 
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The interesting feature which emerges from this calculation is the observed balance be- 
tween the Goldstone contribution and the Quark Loop contribution. Indeed, as the con- 
stituent quark mass Mq gets larger and larger, the Goldstone contribution dominates; while 
for Mq smaller and smaller it is the Quark Loop contribution which dominates. This is il- 
lustrated by the plot of the total aj."^^^^^ (CxQM) versus Mq shown in Fig. (14). What this 
plot shows is in flagrant contradiction with the results reported in ref. [40] based in a calcula- 
tion using a Dyson-Schwinger inspired model. In this model, the authors find a contribution 
from the vr'^-exchange which, within errors, is compatible with the other phenomenological 
determinations and, in particular, with our CxQM result in Eq. (4.16); yet their result for 
the equivalent contribution to the quark loop turns out to be almost twice as large with a 
total contribution 

^(HLbyL) _ [4Qj^ ^ ^217 ± 91) X 10~" . (4.21) 

The central value of this result would require a ridicously small value of Mq in order to be 
reproduced by the CxQM and, furthermore, for such a small value of Mq the vr'^-exchange 
contribution would be far too small as compared to all the phenomenological estimates, 
including the one in ref. [40]. We conclude that a range of values such as 
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Figure 14. 

^(HLbyL) ^^^Q-|^^ versus Mq in in 10"^-*^ units. 



170 < aj^^'^y^) X lO" < 308 , (4.22) 

allowed by the result quoted in Eq. (4.21), cannot be digested within the CxQM and in our 
opinion this casts serious doubts about the compatibility of the model used in ref. [40] with 
basic QCD features. 
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5 Hadronic Electroweak Contributions. 



These are the contributions to the muon anomaly which appear at the two-loop level in the 
electroweak sector. They are the ones generated by the hadronic 77^' vertex, with one 7 and 
the Z-boson attached to a muon line, as illustrated by the Feynman diagrams in Fig. (15). 





Figure 15. 

Feynman diagrams with the hadronic ^^Z vertex which contributes to the muon anomaly. 

These contributions are particularly interesting because, a priori, they could be enhanced 
by a large log(M^/m^) factor. However, due to the anomaly-free coupling assignments in 
the Standard Model, there is an important cancellation of UV-scales between the lepton and 
the quark contributions within a given family [41, 42]. What is left out of this cancellation 
in the sector of the u, d and s quarks, where the strong interactions play a subtle role at long 
distances, is governed by the dynamics of spontaneous chiral symmetry breaking [41, 43-46]. 
The CxQM, where the hadronic blob in Fig. (15) is replaced by a constituent quark loop as 
illustrated in Fig. (16), offers a simple way to estimate these contributions which we next 
discuss. 

In full generality, the hadronic 77 .Z' contribution to the muon anomly, which we denote 
by is given by the following representation [43] : 



,(HEW) 



lim 



d^ 1 



16 cos2 9w Ml fc2^o j (27r)4 - Ml 



kp + 



pi' 



Ml 



r 



— 2q-p 



7 75 + 7 75 — rr^ — ' 

q'^ + Iq-p 



W^,p{q,k), (5.1) 



where W^yp{q,k) denotes the hadronic Green's function: 

W^,p{q,k)= [ d^xe'^i-^ I dSe'^^~'iyy{Q\T{V-{x)Al%y)V-{mQ) 



(5.2) 
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with k the incoming photon four-momentum associated with the classical external magnetic 
field, and where 

V;"'{x)=q{x)^^Qq{x), and A^^ {y) = q{y)j,-f5Q''^\ q{y) (5.3) 

with 

Q = <3L = Qfl = diag (2/3,-1/3,-1/3), and Q^^^ = diag (1, -1, -1) . (5.4) 

The relevant question here is the contribution to a^j^^^^ from the non-anomalous part of 
W^up{q,k), denoted by W^up{q,k), i.e. 

W^,piq,k) = -i-^^^^^^e^p^^q^k^ + W^,p{q,k) , (5.5) 

where the first term in the r.h.s. is the one generated by the VVA anomaly. The second 
term W^i,p{q, k), in the chiral limit where the light quark masses are neglected, is then fully 
transverse in the axial neutral current (i^ index) and the Ward identities constrain it to have 
the form (Q^ = _g2) [43]. 

Wf,up{q, k) = ik" [qpef.yaaq"' - q^epuapq"] W{Q'^) , (5.6) 
with only one invariant function W^Q"^) which depends on the details of the dynamics. 

X 



+ Perm. 




Figure 16. 

Feynman diagrams in the CxQM of the 77.^ vertex type. 
In the CxQM, with = 1, the function W{Q'^) is given by the expression: 

^xQM(g ) = Y^3^ X dxj^ dy-^, 



Q Jo Jo l + ^y{l-y) 



-iVc2 1 
127r2 3 Q2 



l + T#^^log V ^ y (5.7) 

1 + ^ + 1 
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Not surprisingly, the second term in the brackets coincides with the analytic expression of the 
CxQM vertex function given by the term in brackets in the second line of Eq. (4.10). This is 
because, up to an overall factor, W^QyiiQ"^) is precisely the same function as -F^o?^^^* {q'^j Q^j O) 
with the anomaly {Mq — )• oo) subtracted. It is this fact that guarantees that W^qm{Q'^), for 
= 1, has the correct pQCD leading short-distance behaviour [44, 46]: 

W.QmW^) ^ ^1^. (5.8) 

The CxQM, however, has its limitations and does not predict correctly the subleading short- 
distance behaviour: 



Q 



— N 2 



1 Q'n(^^Q^ 



(5.9) 



which falls as O while the OPE in QCD predicts that this subleading term must fall 

as O (^) [43]. 

Another interesting limit is the long-distance behaviour of the function W{Q'^) which in 
QCD is related to a coupling constant of 0{p^) in the odd-parity sector of the effective chiral 
Lagrangian [46] . Here, the prediction of the CxQM is 

W,«m(Q') ^ (5.10) 

Unfortunately, there is no model independent prediction for 1^(0) to compare with. 

The contribution to a!p^^^^ from the anomalous term in Eq. (5.5), evaluated in the 
Feynman gauge [41, 43] is: 



(HEW) 



G^mJaiV, 4 M| 2 
V2 87r2^ 3 ] 3 ^ m2 ^ 3 



ml Ml 




^1 a 

87r2 TT 



X 18.69, 



(5.11) 



and the one from the transverse component W^j_yp{q,k) in Eq. (5.5), evaluated in the CxQM 
in the Feynman gauge and with gA = 



a 



(HEW) 



xQM _ Gf a 

transv 

^/2 87r2 vr 3 



2 Ml 



+ o 



Mi 



Ml 



(5.12) 



The sum of these two contributions takes care of the hadronic sector induced by the dynamics 
of the light quarks u, d and s; but, as already mentioned, it is only when added to the lepton 
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contributions from the electron and the muon and the one from the heavy charm quark that 
the whole sum is gauge independent and it then makes sense in the Standard Model. We 
reproduce below the details of this overall result: 



a 



(HEW) 



-31og^ - - + 21og^ + 1+ log-Tf-2 



mf 



mf 



electron 



u,d anom 



u,d transv CxQM 



-3 log 



Ml 



11 8 

h -vr 

6 9 



2 
3 



log 



Ml 1 1. M2 
2" "3 ~3 



m 



logT^ + o +41og^^ . 5.13 
MA 3 



s anom 



s transv CxQM 



charm 



Notice how the logM| dependence cancels in each generation [45], and we finally obtain 



(HEW) 



y/2 Svr^ TT I 3 



M, 



mf 



A 1 

41og^ 



5 + -vr^ 



a/2 Svr^ TT 



^ - X 18.6 



-5.0 X 10 



-11 



(5.14) 
(5.15) 



for rric = 1.5 GeV and Mq = 240 MeV, a result which within the systematic errors of the 
model, is compatible with the phenomenological determination [45]: 



a 



(HEW) 



e,u,d;ii,s,c 



(-6.7 ±0.5) X 10' 



-11 



(5.16) 



6 Summary and Conclusions. 



From the previous considerations we conclude that the CxQM provides a useful and simple 
reference model to evaluate the hadronic contributions to the anomalous magnetic moment 
of the muon. The effective Lagrangian of this model is renormalizable in the Large-Nc 
limit [12] and, as shown in [11], the number of the required counterterms in this limit is 
minimized for a choice of the axial coupling: = 1. The only free parameter of the model 
is then the mass of the constituent quark mass Mq which in Section II, from a comparison 
with the phenomenological determination of the lowest order hadronic vacuum polarization 
contribution to the muon anomaly, has been fixed to 



Mq = (240 ± 10) MeV , 



(6.1) 
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This range of values for Mq reproduces the phenomenological determination within an error 
of less than 10%. All the other hadronic contributions have then been evaluated for this range 
of values of AIq with the results which are summarized in Table 2. 

Table 2. Summary of results for the hadronic contributions to the muon anomly in the CxQM. 



Class 


Result in 10 units 


HVP 


+47 

652 

-42 


HVP to O (^)^ 


-6.4 ± 1.2 


HLbyL 


15.0 ±0.3 


HEW 


-0.5 



We want to emphasize that the errors quoted in Table 2 are only those generated by the error 
of Mq in Eq. (2.15) and they do not reflect the systematic error of the model. These results, 
within a systematic error of 20% to 30%, are in good agreement with the phenomenological 
determinations. One exception, discussed in detail in Section III.4.1, is the contribution from 
the TT^-y intermediate state to hadronic vacuum polarization where the CxQM, because of its 
SU{3) invariance, fails to reproduce the phenomenological determination which is particularly 
enhanced because of the large observed branching ratio in Eq. (3.54). 

Ironically, the error ±0.3 for the Hadronic Light-by-Light contribution appears to be 
the smallest relative error. This is due to the fact that, as shown in Fig. 16, the sum of 
the quark loop contribution and the Goldstone exchange contribution for values of Mq in 
the range of Eq. (2.15) is already very near to the minimum in the Mg-dependence of the 
sum of these two contributions, which occurs at Mq ~ 300 MeV. In other words, in the 
CxQM the contribution to the muon anomaly which is less sensitive to the value of the 
constituent quark mass is precisely the one from the hadronic light-by-light scattering. This 
fact, however, should not mask the intrinsic systematic error which has not been included. 
Within an expected systematic error of ~ 20%, our results agree with the phenomenological 
determinations reviewed in ref. [39]. An exception, however, is the determination quoted in 
Eq. (4.21). As discussed in the text, the large range of values allowed by this result, cannot 
be digested within the CxQM and in our opinion casts serious doubts about the compatibility 
of the model used in ref. [40] with basic QCD features. 

Another interesting feature, which has appeared when evaluating the hadronic elec- 
troweak contributions, is the impact of the choice = ^, which was initially made on 
theoretical grounds. It turns out that it is only for this choice that the CxQM has the cor- 
rect matching at short-distances with the one predicted by the OPE in QCD [44, 46] when 
evaluating the hadronic electroweak contribution. 

Concerning the next-to-leading contributions from Hadronic Vacuum Polarization we 
have made two observations which are in fact model independent. On the one hand we have 
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explained why this contribution is smaller than the naive expected order of magnitude and on 
the other hand we have derived a sum rule in Eq. (3.27) which offers an interesting constraint 
when evaluating radiative corrections. 
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